Performance of Jacobi preconditioning in Krylov subspace solution of finite element equations     

F.‐H. Lee  K.K. Phoon  K.C. Lim  S.H. Chan 《国际地质力学数值与分析法杂志》2002,26(4):341-372

This paper examines the performance of the Jacobi preconditioner when used with two Krylov subspace iterative methods. The number of iterations needed for convergence was shown to be different for drained, undrained and consolidation problems, even for similar condition number. The differences were due to differences in the eigenvalue distribution, which cannot be completely described by the condition number alone. For drained problems involving large stiffness ratios between different material zones, ill‐conditioning is caused by these large stiffness ratios. Since Jacobi preconditioning operates on degrees‐of‐freedom, it effectively homogenizes the different spatial sub‐domains. The undrained problem, modelled as a nearly incompressible problem, is much more resistant to Jacobi preconditioning, because its ill‐conditioning arises from the large stiffness ratios between volumetric and distortional deformational modes, many of which involve the similar spatial domains or sub‐domains. The consolidation problem has two sets of degrees‐of‐freedom, namely displacement and pore pressure. Some of the eigenvalues are displacement dominated whereas others are excess pore pressure dominated. Jacobi preconditioning compresses the displacement‐dominated eigenvalues in a similar manner as the drained problem, but pore‐pressure‐dominated eigenvalues are often over‐scaled. Convergence can be accelerated if this over‐scaling is recognized and corrected for. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

电离层总电子含量时空特征分析及分区建模   总被引：1，自引：0，他引：1  

刘志平  赵自强  郭广礼 《武汉大学学报(信息科学版)》2012,37(11)

利用地统计方法研究了全球VTEC时空分区特征,提出了IONEX单天产品发布的改进策略和VTEC时空分区建模方法,包括组合函数与IRI2007模型背景场建模方法。纬度10°～70°、经度60°～180°范围内单天VTEC分区建模与结果分析表明,分区建模极大地降低了模型的复杂性。上述两种方法均可得中低纬分区B13约为高纬分区A23中误差（2～3TECU）的2倍,验证了顾及VTEC时空特征分区建模的有效性。  相似文献   

基于QTM的海平面上升分析与模拟   总被引：1，自引：0，他引：1  

侯妙乐  姜晓轶  赵学胜  邢华桥 《地理与地理信息科学》2012,28(1):35-38,2

针对海平面上升影响范围分析与模拟这一国际前沿问题,通过研究基于球面四元三角网(QTM)的关键技术问题,包括复杂拓扑关系计算、LOD剖分、球面水淹分析、基于QTM的多分辨率的DEM数据组织方法和分析精度的相关性评定等,以.Net和Direct3D为开发工具,设计开发了基于QTM的海平面上升影响范围评估模型。该研究结果可为全球海平面上升影响的防灾减灾决策提供有效支撑,并推动了球面数据模型和球面格网拓扑分析的理论成果在全球变化预测相关领域的应用进展。  相似文献   

Temporal and Spatial Variations of the Aerodynamic Roughness Length in the Ablation Zone of the Greenland Ice Sheet   总被引：2，自引：1，他引：1  

C. J. P. P. Smeets  M. R. van den Broeke 《Boundary-Layer Meteorology》2008,128(3):315-338

To understand the response of the Greenland ice sheet to climate change the so-called ablation zone is of particular importance, since it accommodates the yearly net surface ice loss. In numerical models and for data analysis, the bulk aerodynamic method is often used to calculate the turbulent surface fluxes, for which the aerodynamic roughness length (z ₀) is a key parameter. We present, for the first time, spatial and temporal variations of z ₀ in the ablation area of the Greenland ice sheet using year-round data from three automatic weather stations and one eddy-correlation mast. The temporal variation of z ₀ is found to be very high in the lower ablation area (factor 500) with, at the end of the summer melt, a maximum in spatial variation for the whole ablation area of a factor 1000. The variation in time matches the onset of the accumulation and ablation season as recovered by sonic height rangers. During winter, snow accumulation and redistribution by snow drift lead to a uniform value of z ₀≈ 10⁻⁴ m throughout the ablation area. At the beginning of summer, snow melt uncovers ice hummocks and z ₀ quickly increases well above 10⁻² m in the lower ablation area. At the end of summer melt, hummocky ice dominates the surface with z ₀ > 5  ×  10⁻³ m up to 60 km from the ice edge. At the same time, the area close to the equilibrium line (about 90 km from the ice edge) remains very smooth with z ₀ = 10⁻⁵ m. At the beginning of winter, we observed that single snow events have the potential to lower z ₀ for a very rough ice surface by a factor of 20 to 50. The total surface drag of the abundant small-scale ice hummocks apparently dominates over the less frequent large domes and deep gullies. The latter results are verified by studying the individual drag contributions of hummocks and domes with a drag partition model.  相似文献   

Fully coupled hydraulic fracture simulation using the improved element partition method     

Shujun Peng  Zhennan Zhang  Jianye Mou  Bing Zhao  Zhiyuan Liu  Jiading Wang 《国际地质力学数值与分析法杂志》2019,43(1):441-460

The improved element partition method (IEPM) is a newly developed fracture simulation approach. IEPM allows a fracture to run across an element without introducing extra degrees of freedom. It can also simulate any number of fractures in a prescribed mesh without remeshing. In this study, the IEPM is extended to hydraulic fracture simulation. First, the seepage and volumetric storage matrix of a cracked element are derived using virtual nodes (the intersection points of a crack with element edges). Subsequently, the fully coupled hydromechanical equation is derived for this cracked element. To eliminate the extra degrees of freedom (virtual nodal quantities), the water pressure and displacement of the virtual nodes are associated with their adjacent nodes through least squares interpolation. Finally, the fully coupled equation in terms of nodal quantities is obtained. The verification cases validate the method. By using this method, the field-scale hydraulic fracturing process is well simulated. The proposed approach is simple and efficient for field-scale hydraulic fracture simulation.  相似文献   

基于分区间泛化Kriging近似模型的船舶局部结构形状优化     

张干锋  王德禹 《海洋工程》2018,36(6):46-57

为提高Kriging近似模型在船舶结构性能多维度响应预测方面的适用性,对常规Kriging近似模型进行分区间泛化改进:一是对设计样本点的各个维度(分量)进行划区,并在每一个分区间内采用最优拉丁超立方(OLhd)取样;二是引入比例系数w1组合高斯与指数型相关函数提高模型对数据的泛化能力,从而在每个划分的区间内建立泛化的Kriging近似模型。通过SCH测试函数,验证了构造的分区间泛化Kriging模型有效性。结合参数化建模和改进Kriging近似模型对某过渡肘板和舱口角隅边界进行形状优化,优化第一阶段由改进的近似模型通过多岛遗传算法得到全局初步的最优解,第二阶段在初步解的基础上缩小优化变量范围,由少量的FEM计算即可在小范围内搜寻到精确的最优设计变量。结果表明:分区间泛化Kriging近似模型在预测多维度响应时较常规Kriging模型预测精度更高;分阶段的形状优化流程在保证极小误差和缩小计算成本的情况下能够得到理想的应力分布和重量优化结果,有助于船舶结构的轻量化研究。  相似文献   

Quantifying sand delivery to deep water during changing sea-level: Numerical models from the Quaternary Brazos Icehouse continental margin     

Oriol Falivene  Bradford E. Prather  John Martin 《Basin Research》2020,32(6):1711-1733

Sequence stratigraphy for clastic continental margins predicts the development of sand-rich turbidite deposits during specific times in relation to base-level cycles. It is now widely understood that deltas can extend to the shelf-edge forced by high sediment flux and/or base level, providing a direct connection to transfer sediment and sand to the slope and basin floor even during high base level periods. Herein, we build a stratigraphic forward model for the last 120 kyr of the fluvio-deltaic to deep-water Brazos system (USA) where sediment partitioning along an Icehouse continental margin can be evaluated. The reduced-complexity stratigraphic forward model employs geologically constrained input parameters and mass balance. The modelled architecture is consistent with the location of depositional units previously mapped in the shelf. Sand bypasses the shelf and upper slope between 35 to 15 kyr before present and only about 20%–30% of all the sediment and sand supplied to the system is transferred to deep water. Several scenarios based on the initial Brazos model investigate the relationships between base level and deep-water sand ratio (DWSR). DWSR is defined as the relative amount of sand transferred to the deep-water portions of the system subdivided by the total sand input to the model. Linear correlations between DWSR and base level change rates or base level are very poor. Short-term variability due to local processes (for example avulsions) is superimposed to the long-term trends and mask the base level signal. DWSR for an entire base-level cycle is mainly controlled by the proportion of time the delta stays docked at the shelf-edge. Stratigraphic forward models are useful to complement field observations and quantify how different processes control stratigraphy, which is important for making predictions in areas with limited information.  相似文献   

Isolation of plane shear wave using water saturated trench barrier     

《Soil Dynamics and Earthquake Engineering》2014

The article presents a theoretical investigation of the propagation of normally and obliquely incident plane shear waves past a rectangular trench filled up with water saturated soil sandwiched between anisotropic elastic media. The motivation for this work is due to the effectiveness of the isolation of waves by the in-filled rectangular trench. Unlike the most of the previous researchers, this model considers the soil–structure interaction effects and directly determines the influence of barrier in the form of in-filled trench on the mode of wave propagation. It is of interest to determine the reflection and transmission coefficients, and the energy partition distribution of shear waves in the in-filled rectangular trench showing the influence of barrier on the propagation of waves. An extensive parametric study through numerical computation is carried out to investigate the influence of the material properties of the in-filled trench and the amplitude ratios on shear waves. The in-filled trench barrier directly declines the intensity of waves significantly in such a way that the waves do not create any hazards to the nearby structures, if exists at all.  相似文献   

Structural modal parameter identification and damage diagnosis based on Hilbert-Huang transform   总被引：1，自引：1，他引：0  

Han Jianping  Zheng Peijuan  Wang Hongtao 《地震工程与工程振动(英文版)》2014,13(1):101-111

Traditional modal parameter identifi cation methods have many disadvantages,especially when used for processing nonlinear and non-stationary signals.In addition,they are usually not able to accurately identify the damping ratio and damage.In this study,methods based on the Hilbert-Huang transform(HHT) are investigated for structural modal parameter identifi cation and damage diagnosis.First,mirror extension and prediction via a radial basis function(RBF) neural network are used to restrain the troublesome end-effect issue in empirical mode decomposition(EMD),which is a crucial part of HHT.Then,the approaches based on HHT combined with other techniques,such as the random decrement technique(RDT),natural excitation technique(NExT) and stochastic subspace identifi cation(SSI),are proposed to identify modal parameters of structures.Furthermore,a damage diagnosis method based on the HHT is also proposed.Time-varying instantaneous frequency and instantaneous energy are used to identify the damage evolution of the structure.The relative amplitude of the Hilbert marginal spectrum is used to identify the damage location of the structure.Finally,acceleration records at gauge points from shaking table testing of a 12-story reinforced concrete frame model are taken to validate the proposed approaches.The results show that the proposed approaches based on HHT for modal parameter identifi cation and damage diagnosis are reliable and practical.  相似文献   

微量元素在湖积物、土壤的垂向分布与稻谷中的分配   总被引：22，自引：7，他引：15       下载免费PDF全文

廖启林  黄顺生  范迪富  吴新民  朱伯万  翁志华  华明  陈宝  冯金顺 《第四纪研究》2005,25(3):331-339

报道了“江苏省国土生态地球化学调查”项目有关湖积物、土壤柱微量元素垂向分布及稻谷样元素分配研究的初步成果。通过对有关样品元素含量的分析对比,指出洪泽湖沉积柱元素垂向分布与其动荡的沉积环境关系密切、太湖局部沉积物在其目前20cm以上深度存在人为重金属污染,这二湖沉积物中其氮与总有机碳含量呈显著正相关;认为城市化可使当地土壤环境受重金属等污染的深度显著加大、城市土壤环境受人为活动的影响要显著高于农村;发现Mo易富集于稻米、重金属与营养元素易富集于稻皮、各元素在米-谷间的分配关系不确定。  相似文献